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ABSTRACT 

In practice, thyristor rectifiers are used as high 
current-low voltage power supplies. Under certain con- 
ditions such rectifiers introduce subharmonic resonances. 

In this paper, a limit cycling voltage regulator 
system is simulated in the digital computer. The ref- 
erence voltage and load are varied and the system gain 
is changed for better understanding of .their effects on 
the subharmonics. A Fourier analysis of the output wave 
shape then determined the existence and amplitude of 
subharmonic components. 
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I. LIMIT CYCLING SYSTEMS, THEORY OF OPERATION 



A. INTRODUCTION TO LIMIT CYCLING SYSTEM 

In practice 3 it might be required to have a control 
system which will generate a limit cycle. In this thesis 
the goal is to investigate the characteristics of a limit 
cycling voltage regulator. This kind of voltage regulator 
system is potentially useful as a low cost source of well 
regulated direct current 

The system introduced in this paper has linear and 
nonlinear elements. One of the tools for linearizing non- 
linear systems is the describing function method. In deal- 
ing with limit cycling systems the ordinary describing 
function method predicts the limit cycle a but is not ade- 
quate for a complete analysis . The Dual input describing 
function is a useful method for investigation of the limit 
cycling system when subharmonic resonances are anticipated. 

B. DUAL INPUT DESCRIBING FUNCTION (DIDF) 

West (^) 3 introduced a new mathematical device which 
is called the "Dual input Describing function" (DIDF) . It 
is an extension of the conventional describing function 
which linearizes the nonlinearity when forced by two inputs. 

The Dual input Describing function is a very powerful 
tool for the investigation of limit cycling systems. It 
will be used for prediction of subharmonic resonance. 
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For better understanding of a limit cycling system a 
general block diagram is introduced in Figure 1-la. In 
this limit cycling control system N(A*w) represents the non- 
linear element, G(ja)) the linear element. The DIDF is 
based on the error input to the nonlinear element. For a 
well designed system the error X(t)* is equal to 

X(t) = r(t)-C(t) 

and would be small so that any limit cycling period To* can 
be modeled as shown on Figure 1-lb. In this model the sin- 
usoid amplitude A is associated with. one input and the D.C 
bias amplitude is associated with the second input. Then 
the input to the nonlinear element is 

X(t) = B+Asine 

Now* the calculation of the DIDF for the nonlinear element 
N(a*co)* is accomplished by computing the ratio of the out- 
put amplitude spectrum of the nonlinearity to its input 
which gives two DIDF. 

N (A 3 B) = e” J0 
° A 

and 

N s ( A, B) = 

The first nonlinearity is considered as limit cycle DIDF 
and will be used to predict the system limit cycle. The 
second DIDF is considered as a linear gain for range B/A 
and leads to linear system description due to input. The 
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FIGURE 1.1a. Basic limit cycling system. 




FIGURE 1.1b. Representation of input-output relation of a 
nonlinear element. 
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OUTPUT 




FIGURE 1.1c. Input and output wave form for an ideal relay. 
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mathematical definition of both DIDF’s follow directly from 
Fourier analysis and they can be written as 

N o (A,B) = Re[N 0 (A,B)]+JI m [N 0 (A,B)] 



1 

ttA 



2 7T 

J y (B+Asin0 )sin0d0 
o 



+ 



J 

ttA 



; 



y (B+Asin0)cos0d0 



and 



N s (A,B) 



1 

2ttB 



2 tt 

J y (B+Asin0 ) d0 
o 



( 1 - 1 ) 



(1-2) 



Equation (1-1) can be reduced since integration of the co- 
sine term of N q (A 3 B) in this interval is zero. Then the 
simplified form of N (A,B) can be written as 



n o (a 3 b) 



1 

2ttB 



2 7T 

J y (B+Asin0 ) sin0d0 
o 



( 1 - 3 ) 



In the limit cycling system B/A will be small .and an incre- 
mental DIDF can be derived as B->o 



DIDF A lim [N (A,B)] 
B+o S 



and 



DIDF 



1 lim 

2tt B->-o 



_ 1 _ 

B 



J 



y (B+Asin0 )d0 



( 1 - 4 ) 
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C. DIDF FOR IDEAL RELAY 

Let us consider an ideal relay with drive levels ±D, 
Figure 1.1c. There the input to the nonlinearity is 

X = B+Asino) 0 t 
= B+Asin0 



The ideal relay has a single valued nonlinearity, there- 
fore, computation of the limit cycle DIDF is done by using 
the result obtained in equation (1-3) and applying it to 
Figure 1.1c. Then the DIDF can be written 



and 



27T 

q (A,B) = ■ — - J y (B+Asin0 ) sin0d0 



7T+01 27T-0! 

j- J Dsinede + J (-D)sinede 



ttA 



7T + 0 ] 



2 TT 

+ J Dsin0d0 

2 TT — 0 i 



4D 

ttA 



cos 0 1 



6i 

N 0 (A,B) 



_1 

sin (B/A) 

4P , 

ttA 1 - (B/A) 2 



( 1 - 5 ) 
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The DIDF 3 associated with the signal is 



n s (a 3 b) 



1 

2ttB 




Dd0 



2TT-0 




2TT 

(-D)d + J Dd0 

27T-0! 



2D 

ttB 



01 



= -^g- sin 1 (B/A) 

If B+o , forumla (1-6) can be written 

N (A,B) = lim [N (A ,B) ] 
s B+o s 



( 1 - 6 ) 



N g (A,B) = 



2D 



irA (1-7) 

And the DIDF for two sinusoids input is found as: 



X = Asin(a) t + 0)+Bsinto t 
o s 



in which o) q is the limit cycle frequency and o) s is the sig- 
nal frequency. Then the DIDF will be 



N g (A,B) 



N s (A,B) 







, X 


2 
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4D 


n 1 
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B 




ttA 




A, 


" CT 


Aj 
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2D 


1 1 




I 2 x 3 1 




|\ 


IT 


1+ 8 1 


UJ 


1 + CT | 
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| +... 



( 1 - 8 ) 



(1-9) 



If this result is compared with the approximate DIDF, it can 
be said that the approximate DIDF which is derived as equa- 
tion (1-5), (1-6) are valid for 5 percent accuracy but if 
the following conditions are satisfied: 
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( 1 - 10 ) 




Then it is possible to have 2.5 percent accuracy with re- 
spect to the exact DIDF. 

Therefore, for the ideal relay case the limit cycle 
DIDF are always valid. to better than 5 percent accuracy re- 
ferred to the exact DIDF, under the restriction of (1-10). 

D. LIMIT CYCLING SYSTEM CHARACTERIZATION 

The system here considered is a kind of system which 
always has a limit cycle. Most limit cycling systems in- 
clude a pre-filter like that shown in Figure 1.2. The 
fundamental purpose of this filter is to reshape the command 
signal amplitude to the input such that at the input of the 
nonlinearity inequalities (1-10) are satisfied. This brings 
a design problem of the pre-filter which will not be con- 
sidered in this thesis. 

1 . Frequency Domain Description 

To investigate sinusoidal frequency analysis of 
limit cycling systems, the system model is illustrated in 
Figure 1.2. r(t) is a sinusoid and the input to the non- 
linearity N is 

X = Asin(u) Q t+0 )+Bsina) s t 

The goal here is to investigate the interaction between two 
sinusoids, through the nonlinearity and differents inputs 
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are considered in this analysis. Figure 1.3 shows the range 
of possibilities of the signal and limit cycle. For Fig- 
ure 1.3a the output through the nonlinearity is nearly sin- 
usoidal at the limit cycle frequency but is phase shifted 
with respect to the limit cycle. It can be said that sig- 
nal through relay is zero which means N g =0. Also it may 
be said that for signal frequency greater than the limit 
cycle frequency but satisfying the amplitude ratio con- 
dition also gives N s =0. 

In Figure 1.3b the amplitude ratio is satisfied 
but the frequency ratio is 0.5* In this case the signal 
through the nonlinearity would not lose identity and the 
amplitude ratio only will depend on the phase shift of the 
input component to the nonlinearity. Then for the case 
shown in Figure 1.3b, it can be concluded that N s =N s (0). 

For a frequency ratio equal to 0.33, affects the input 
signal amplitude by factor of 16 percent and a phase of 
4 degrees. For a frequency ratio greater than 0.33, (but 
amplitude ratio satisfied) N g can be considered independent 
of relative time position since the change in amplitude and 
phase ratio are never greater than 1 percent. Therefore as 
a result of these conclusions the open loop transfer function 
can be written as 

co 

G eq ( ' 3lj)=N s- H eq- G(J “ ) for |u|< T 

CO 

G e q(J“) = T-N s . H eq- G (J“) for j~ < M <u 0 

?eq (Ja))=0 for Vl“l (1-11) 
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FIGURE 1.3a. AR = 0.33, FR = 1.00 





FIGURE 1.3. Range of possibilities of input. 
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where y represents a complex number depending on the rela- 
tive phase of these signals . 

2 . Pole-Zero System Characterization 

Frequency domain analysis shows that in the range 

between -a) o /3<oo s <u) o /3 rad/sec along the imaginary axis, if 

the amplitude ratios are also satisfied the system behaves 

linearly. In the ranges to /3<w <w and -to <to <to /3 rad/sec 

o - s- o o- s- o 

the signal gain of the nonlinearity is a function of time, 
and for the remaining portion of the jco axis, the signal 
gain is always zero. Now, consider exponentials for which 
the amplitude ratio is satisfied, (amplitude ratio refers 
to ratio of peak exponential to peak limit cycle amplitude) 
it may be noted that the exponential excited output of a 
relay is relatively independent of the time origin of the 
input exponential signal in the existence of the limit cy- 
cle, if the signal lasts for two or more limit cycle periods 
Figure 1.4 shows this situation. Then, the minimum value 
of the input time constant is, 



Then, the pole-zero configuration of the system may be 
described. Consider the complex plane "p" to be divided 
into three regions, shown in Figure 1.5. 

The nonlinear element N of the system will be 
considered as a gain element and N is derived for the relay, 
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FIGURE 1.5. Root locus plot for G(p) . 
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therefore, it will not contain any phase shifting. If the 
open loop system gain is adjusted to one-half of the value 
of the limit cycle gain, it is possible to automatically ob- 
tain a limit cycling system. The closed loop roots of the 
system can be approximately determined, such that, if the 
gain at the crossover of the imaginary axis N=N q which is 
the limit cycle gain, and at gain value of N=N q /2=N s one 
of the closed loop roots will be close to the intercept 
point of the w q /3 circle on the real axis which is shown as 
P & . anc ^ k* 16 °^ er complex pair of roots would be 

very close to the circumference of the co circle. This 

o 

result determines the dynamic response of the limit cycling 
system . 

One of the important points about limit cycling 
systems is from the stability point of view: if the sensi- 

tivity of the system is increased then the limit cycle am- 
plitude increases but the gain N g decreases and always the 
S.N g product remains constant. Closed loop root locations 
do not change, only the limit cycle amplitude will be 
changed. Also, it is possible to have instability in the 
system but may be the system response does not change. 

E. CLOSED LOOP PROPERTY 

Consider a limit cycling system as shown in Figure 1.2. 
Then the open loop transfer function of the system is 

G eq (p) = N s -H 1 (p)-H 2 (p).G(p) 

= N s .H eq .G(p) (1-12) 
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The DIDF for the ideal relay were found in equations (1-5) 
and (1-6) and they may be approximated as, 

N * 

O IT A 



N « 
s ttA 



and the ratio of these two. 



N 

s - 1 

N 2 

• o 



then N g may be written in equation (1-12) 



G (jco) = N .H .G(jaj) 
eq J s eq w 



1 N .H . G ( j to ) 
= o ° ec l 



in order to sustain the limit cycle 



N 0 - H eq- G (J“) " - 1 



(1-13) 



must be satisfied. Then, 






( 1 - 1 ^) 



This is a kind of proof of the explanation given earlier 
for open loop gain of limit cycling systems. 
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II. PERFORMANCE CHARACTERISTICS OF THE LIMIT 
CYCLING VOLTAGE REGULATOR 



A. SYSTEM DESCRIPTION 

The block diagram of the voltage regulator system used 
in this thesis is shown in Figure 2.1. The complete circuit 
diagram is shown in Appendix A. 

The firing circuit, instead of using a separate cir- 
cuit for each thyristor uses a relay to turn all thyristors 
on or off. If the relay is on a positive voltage is applied 
to all thyristor gates and this causes conduction in the 
rectifier. If the relay is off then the rectifiers will 
not conduct. . 

The rectifier unit is a three phase full-wave recti- 
fier made of a combination of six thyristors. The thyristor 
is a solid state device which is equivalent to a thyratron. 
One thyristor device is a silicon controlled rectifier. It 
is a four layer device and may be likened to a PNP transis- 
tor joined to an NPN transistor and sharing- a common col- 
lector-base junction. The thyristor and the characteristic 
curve of the thyristor are shown in Figure 2.2. 

As a load unit, a second order filter with very low 
damping factor is used. It is an LCR filter. The values 
of L,C are constant and variation of the damping ratio is 
obtained by varying the load resistance. 

The system has a nonlinearity due to the relay and 
the thyristor characteristics provide a second nonlinearity. 
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FIGURE 2.1. Basic block diagram of the limit cycling 
voltage regulator. 
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FI.GURE 2.2. Thyristor characteristic. 
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These two nonlinearities are treated as one nonlinearity in 
analysis of the system^ providing the simplified block dia- 
gram of the system as shown in Figure 2.3. 

B. STEADY STATE ACCURACY 

It is important to investigate the system response to 
different inputs. For this purpose the voltage regulator 
system is analyzed for a load resistance R=0.8ft, and dif- 
ferent reference voltages applied to the system. The im- 
portant feature to be studied is the output characteristic 
of the system. For different reference voltages the out- 
puts obtained were subjected to Fourier analysis. Table 2-1 
shows the result of this analysis. In Table 2-1 results 
taken from digital computer calculation of the Fourier se- 
ries and subroutine FORIT used for this analysis. In the 
Table 2-1 A q is the DC value of the output. 

These results show that the system output is not equal 
to the required value. Then the difference may be defined 
as a a 

a = V « - A 
ref o 

For high reference voltages the system has large a and if 
Vref gradually decreases 3 a also decreases. At ^ re f = 52 
volt a changes it's sign 3 that is, for low voltage cases 
the system output voltage is a little bit higher than V re f 
This may be good but when decreases the system ripple 

instability increases which is not desired for the system. 

The numerical results show that if ^ re f is decreased gradually, 
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Simplified block diagram of the limit cycling 
voltage regulator. 



FOURIER ANALYSIS FOR R=0.8n 



voltage 

VrefM 


A„ Ivolt) 


LIMIT CYCLE 
AMPLITUDE 
(volt) 


SUBHARfWiC 

Amplitude 

(Volt) 


LiMiT CYCLE 
mQUENCY 

(Hz) 
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ORDER OF 
iUBHARMOWit 
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... 
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... 
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1/2 
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56.277 
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51.2 
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1/2 
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1.53 
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72.4 


36.2 


1/2 


52 


52.966 


2.21 


0.133 


60.4 


30.2 


1/2 


50 


50.272 


1.24 


0.143 


90.0 


45.0 


1/2 


48 


49.181 


1.505 


0.475 


56.2 


28.1 


1/2 


46 


47.497 


1.41 


0.635 


68.4 


17.1 


1/4 



TABLE 2-1 
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subharmonics become more effective and this may make the sys 
tern unstable. For V ^ = 58 volt, the amplitude of the sub- 
harmonic is quite small (about 0.146 volt) but for V ^ = 

46 volt, the subharmonic amplitude becomes 0.635 volt which 
is quite large compared with 0.146 volt. Also the results 
show that the limit cycle amplitude shows an opposite re- 
action. Increase in the subharmonics amplitudes caused a 
decrease in limit cycle amplitude. Figure 2.4 shows the 
output of the voltage regulator at = 60 volt and 

V ref = 58 volt ‘ 

C. TRANSIENT RESPONSE OF THE SYSTEM 

For investigation of the transient response of the sys 
tern only the linear parts are taken into account. The be- 
havior of the nonlinear part is not considered. The linear 
filter of the system is assumed to be as shown in Figure 2.5 
And the transfer function of the system is 



In general second order system characteristic equations can 
be written in the form 

S 2 + 2£co S+o) 2 = 0 
n n 

then for the filter. 




1/LC 



S*+(1/RC)S+(1/LC) (2-1) 



1 



and 




(A). 



n / LC 
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FIGURE 2.4b. Output waveform: v re f = 58V., R = 0 . 8Q 
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FIGURE 2.5. Filter and load of the system. 
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During analysis of the system L and C are held constant and 
only the value of R is changed 



C = O.lf. 

L = 0.0005h. 



therefore , 



a) - 141 rad/sec 
n 

and for R = 0 . 8ft 



£ = 0 .04419*1 

which is a very low damping factor. 

The settling time of the filter is 




Then, 



t 

s 


= 2.5 


sec 


for 


R = 3 . 2 fi 


t 

s 


= 1.23 


sec 


for 


R = i. 6 n 


t 

s 


= 0.645 


sec 


for 


R = 0.8« 



The filter transfer function can be written in Bode form. 



G ( j oj ) 



1 

LC( ju>) 2 +L/R( ja))+l 



And for R = 0.8ft 

G(ja,) = (1/2 10 4 ) ( j(u) 2 + ( 12 . 3/2 10 4 (ju) + l 

(2-2) 

The phase and amplitude curves of the Bode diagram are 
drawn with the IBM-360 digital computer and the curves are 
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illustrated in Figure 2.6. The analysis of the closed loop 
characteristic of the filter is shown on the Nichols plot 
in Figure 2.7* 

D. PREDICTION OF LIMIT CYCLES 

For prediction of the limit cycle frequency, the basic 

describing function method is used. The describing function 

for the nonlinearity of the system which includes relay and 

rectifier is derived by Leszcynski (Ref. 15). And for the 

linear part of the system for this prediction equation -(2-2) 

is used. Figure 2.8 and Figure 2.9 shows the Nichols plot 

of the limit cycling voltage regulator for = 64 volt 

and V ref . = 62 volt both at R = 0.8ft. Figure 2.8 gives the 

limit cycle frequency f^ c = 35HZ and limit cycle amplitude 

B = 3.6544 volt. This agrees with the experimental result 

shown in Table 2-1. Also for V ~ = 62 volt the Nichols 

ref 

chart of Figure 2.9 gave the limit cycle frequency f lc =37HZ 
and limit cycle amplitude about B = 2.4632 volt. Experi- 
mental results in Table 2-1 gives f^ c = 36.4HZ and limit 
cycle amplitude about B = 2.82 volt which agree fairly well 
with the result obtained by Nichols plot. 

E. SUBHARMONICS 

One of the purposes of this paper is also to investi- 
gate subharmonic resonances in the system. Table 2-1 for 
Fourier analysis shows relationships between ^ re f and the 
subharmonic amplitude and frequency. Results show that 
changes in V ^ affect subharmonic amplitudes. If ^ re f is 
increased the subharmonic amplitudes also increase and for 
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FIGURE 2.7. Nichols plot of the filter for R = 0.8S2. 
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FIGURE 2.8. Nichols plot of system for V re£ = 64V., R-0.8B, f lc 35. 
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FIGURE 2.9. Nichols plot of the system for V f =62V. , R=0.8G, 
f lc = 37 HZ. re 
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Vpgf < 50 volt subharmonics become very important. System 
behavior gives a ripple instability at the output. This in- 
stability may be caused by the high amplitude subharmonics. 
Figure 2.10 shows this kind of output for = 48 and 

V ref = 46 volt - 

From the frequency response point of view, it may be 
concluded that changes in V re f do affect appreciably the 
frequency of the subharmonics . The basic frequency of the 
subharmonic tends to remain at 25 cps . Table 2-1 shows that 
changes in the system load resistance; R = 0.8ft, R = 1.6ft, 

R = 3»2ft caused no change in the subharmonic frequency and 
in like manner changes in the reference voltage did not 
change the subharmonic frequency. These results shown in 
Table 5-1, Table 5-2, and Table 5-3* It may be concluded 
that any change in the system will not effect the subhar- 
monic frequency since the system forced the subharmonics 
toward the 25 HZ value. But the amplitude of the sub- 
harmonic is quite sensitive to any change in the system or 
any disturbance to the system. 
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FIGURE 2.10b. Output waveform for R = 0.8ft, V £ - 48V 



III. INTRODUCTION TO THEORY OF SUBHARMONIC RESONSNCE 



Subharmonic resonance is a phenomenon which occurs in 
nonlinear systems. Usually the nonlinear system subharmonic 
resonance frequency is a fraction of the input frequency. 

As an example consider the limit cycling voltage regulator 
system which is introduced in this thesis, or a relaxation 
oscillator, or a moving-coil loudspeaker with a non-uniform 
magnetic field are the common examples of the subharmonic 
resonance. This phenomenon can also occur in servo mecha- 
nisms with saturation or with other nonlinearities. In all 
of these cases the lowest frequency present at the output 
is a submultiple of the input frequency and usually at or 
near the natural frequency of the system. The system out- 
put amplitude would be different depending on the existence 
of the subharmonic. Such subharmonics can increase or de- 
crease the output amplitude. For systems in which the out- 
put is intended to reproduce the input as closely as 
possible, as is the motion of the loudspeaker cone or a 
servo motor, this phenomenon can introduce an intolerable 
distortion. Also it is observed that in limit cycling vol- 
tage regulator systems, subharmonics can cause instability 
of the system. Also one of the important points, subhar- 
monic resonance does not occur in pure linear systems. If 
subharmonics are not wanted the best way to eliminate them 
is by changing the parameters of the linear system or plac- 
ing a compensation element rather than trying to change the 
nonlinear part of the system. 
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It is possible to predict subharmonic resonance in 
the nonlinear system by using the Dual input Describing 
function and by solution of the nonlinear differential equa- 
tion. There are some disadvantages to the second method. 
Since the output of the nonlinear element is expressed as a 
power series of the input, successive approximations to the 
correct solution converge rapidly only if the amount of 
nonlinearity is small. And also the solution of the dif- 
ferential equation does not give a physical picture of the 
mechanism of subharmonic resonance. 

Subharmonic resonance normally exists only in lightly 
damped systems with nonlinear energy delays or with nonlin- 
ear restoring forces. In systems giving rise to subhar- 
monics, the nonlinear element is frequently characterized 
by hysteresis. Typical nonlinearities are backlash, mag- 
netic hysteresis loops and relays, etc. It is also impor- 
tant to note that generation of subharmonics is dependent 
on initial conditions, i.e., they may occur in a certain 
range of amplitudes and frequencies. If in a system a 
subharmonic is obtained it stays "locked in" over a reason- 
able frequency range. If the driving frequency is changed 
to a new value, the system output frequency changes also 
to a new value which is co/n, to being the driving frequency 
and n order of subharmonic. 
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IV. DUAL INPUT DESCRIBING FUNCTION FOR 
PREDICTION OF SUBHARMONICS 

In the voltage regulator system, power supply ripples 
feed around the loop and interact with the nonlinearity of 
the system and cause instability of the system. The limit 
cycle resulting from this instability is at a submultiple 
of the ripple frequency and is commonly called a subhar- 
monic resonance. Normal describing function theory cannot 
predict the existence of the subharmonic resonance. But 
by using the Dual input describing function method, it is 
possible to find a mathematical description of the nonlin- 
ear part of t.he system which can be used to predict the 
existence of these resonances. In general the forward 
transfer function of the voltage regulator system can be 
written 

= N(Ju).G(ju) 

E (4-1) 

and 

N(jco) = N o (jco)+N 1 (jm) 

in which N (jco) is the DIDF due to the DC input to the non- 
linearity of the system and N^(ju) is due to AC input to 
the nonlinearity. G(jm) is the transfer function of the 
linear part of the system. If E represents an input to 
the nonlinearity then equation (4-1) can be rewritten, 

Vout _ N . G(.ico)+N n . G(jo)) 

E = ° (4-2) 
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Frequency analysis is accomplished in the usual fash- 
ion, by applying the Nyquist criterion. The result could 
be interpreted as an encirclement of the -1+jo point. 
Therefore, equation (4-2) can be drawn in vector form as 
shown in Figure 4.1. It is seen from Figure 4.1 that the 
N^.G(joj) vector is placed at the tip of the N Q .G(ju)) vec- 
tor. Because N q ( j co ) represents the DC part of the DIDF 
and N^(jw) represents the AC part of the DIDF, N (ju>) will 
have one component but N^(jw) will have many components as 
the phase angle changes. But during this variation oo = 

must be constant. Then if (for different oj ) this vector 

x 

plot is repeated Figure 4.2 would be obtained. The N^.G(joi) 
vector would be placed on the tip of the N G(jw) as phase 
is varied. By following this procedure many curves can be 
drawn and if any of these curves N^.G(jw) encircles the 
-1+jo point then in Nyquist sense the system may be un- 
stable. Such a system may have a subharmonic resonance, 
but this is not guaranteed. 
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FIGURE 4.1. Vector diagram using DIDF. 
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FIGURE 4.2. Analysis of the DIDF. 



49 



V. ANALYSIS OF THE EXPERIMENTAL RESULTS 



A. ANALYSIS OF LIMIT CYCLING VOLTAGE REGULATOR 

The goal is to find out the exact cause of subharmonics 
in the limit cycling voltage regulator system. For this pur- 
pose the IBM -360 computer was used for solving system dif- 
ferential equations and subroutine INTEG 2 is used. Appen- 
dix B shows this computer program and as an example one of 
the outputs. The output graph was obtained for three dif- 
ferent loads and in each case ^ re f is changed from 64 volt 
to 46 volt. Each run was permitted enough time so that the 
system reached steady state . From these output plots by 
analyzing the wave shape of the output, the results in 
Table 5-1, Table 5-2, and Table 5-3 were obtained. 

Those tables showed that the system has mostly 1/2 
order subharmonics and a few 1/4, 1/8 subharmonics. Then 
it may be said that the limit cycling voltage regulator 
system has an even order nonlinearity. Because most ex- 
perimental results have shown that nonlinearities of even 
order give even order subharmonics. 

In the Tables amplitudes are given peak to peak in 
volts. These amplitudes are given as the observed values 
of the output’s highest spike and lowest spike. These low- 
est and highest spikes become smaller if the load resistance 
is increased. Increasing load resistance caused increase 
in the ripple instability at the output. For many cases it 
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ANALYSIS OF OUTPUTS FOR R=0.8 a 



REFERENCE 
VOLTAGE 
VftEF (VOLT) 


LARGER 

AMPLITUDE 

UP- volt) 


SMALLER 

AMPLITUDE 

(PP.vou) 


LIMIT CYCLE 
FREQUENCY 
(Hi) 


SU3HARMONK 

FREQUENCY 

(Hz) 


. ORDER OF 
SUBHARMONK 


64 


8.66 


8.66 


29.8 


--- 


--- 


62 


7.1 


6.7 


35.8 


17.85 


1/2 


60 


5.9 


5.9 


39.7 


--- 


--- 


58 


4.72 


4.32 


50.8 


25.4 


1/2 


56 


5.5 


5.11 


50.8 


25.4 


1/2 


54 


3.5 


3.15 


68.6 


35.8 


1/2 


52 


4.72 


4.72 


57.8 


--- 


--- 


50 


2.76 


2.36 


87.8 


' 44.6 


1/2 


48 


5.31 


1.97 


56.2 


28.1 


1/2 


46 


11.8 


1.77 


68.4 


17.1 


1/4 



TABLE 5-1 
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ANALYSIS OF OUPUTS FOR R = l.6 fi 



REFERENCE 
VOLTAGE 
VfttF (von) 


LARGER 
AMPLITUDE 
(?L volt) 


SMALLER 

AMPLITUDE 

(?F.volt) 


LIMIT CYCLE 
FREQUENCY 
(Hr) 


SUBHARMONIC 

FREQUENCY 

(Hr) 


ORDER OF 
SUBHARMONIC 


64 


5.3 


5.7 


27.8 


13.35 


1/2 


62 


3.54 


3.54 


35.8 


--- 


... 


60 


1.97 


1.77 


57.8 


28.7 


1/2 


58 


1.18 


1.18 


134.0 


... 


... 


56 


0.985 


0.985 


115.5 


... 


... 


54 


2.76 


0.788 


115.0 


25.4 


1/4 


* 

52 


3.15 


0.984 


145.6 


72.8 


1/2 


50 


3.15 


0.59 


57.8 


25.4 


1/2 


* 

48 


3.55 


1.22 


73.68 


3.21 


1/8 


* 

46 


4.32 


1.37 


80.8 


40.4 


1/2 



* Subharmonic is determined by Fourier analysis. 

TABLE 5-2 
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ANALYSIS OF OUTPUTS FOR R=-3.2 n 



REFERENCE 

VOLTAGE 

VftEf- (VOLT) 


LARGER 

AMPLITUDE 

(P?.Volt) 


SMALLER 

AMPLITUDE 

(pp.volt) 


LIMIT CYCLE 
FREQUENCY 

(Hz) 


SUBHARMONIC 

FREQUENCY 

(Hz) 


ORDER OF 
SU3HARMONI 


64 


3.15 


3.15 


21.2 


--- 


--- 


62 


1.57 


1.48 


39.7 


19.85 


1/2 


60 


1.18 


0.788 


68.8 


36.2 


1/2 


* 

58 


2.36 


0.394 


104.8 


52.4 


1/2 


56 


1.575 


0.985 


90.6 


25.5 


1/2 


* 

54 


1.77 


0.59 


100.0 


50.0 


1/2 


52 


1.575 


0.7 86 


91 


23 


1/4 


* 

50 


5.9 


0.394 


63.6 


17.4 


1/4 


OO 

* 


2.75 


1.18 


62.88 


31.44 


1/2 


* 

46 


3.94 


0.53 


49 . 28 


12.32 


1/4 



*Subharmonic is determined by Fourier analysis. 

TABLE 5-3 
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was not possible to determine the order of the subharmonics 
by inspection of the output curves. In these cases a 
Fourier analysis was made to establish the subharmonic fre- 
quencies and amplitudes . 

Figure 5* la to 5 • Id shows the wave form at the filter 
element and from Figure 5.2 to Figure 5*34 shows the output 
and input plot to the filter. If one of the input plots is 
carefully observed one notices the difference in width of 
these pulses. The system has sinusoidal input to the thy- 
ristors and the output of this rectifier goes to the fil- 
ter. This will charge the capacitor at the filter and when 
the output reaches the relay would be open but at this 

moment the rectifier will still continue to conduct and 
during these times the capacitor will charge. Afterwards 
it will discharge till the output voltage reaches V 
If the output voltage is much higher than V ^ then the ca- 
pacitor will discharge for a longer time, therefore, the 
system would give a different pulse width due to this ef- 
fect. This can be seen by the different number of ripples 
at the top of the input voltage wave form. 

Output plots from Figure 5.2 to Figure 5*11 are for 
load resistance R = 3.2ft. At lower voltages deciding the 
order of the subharmonic by inspection becomes impossible. 
Those outputs were analyzed in Fourier point of view and 
those outputs obtained by using Fourier analysis are desig- 
nated with asterisks in the Table. 
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FIGURE 5.1a.. Input to the filter, V ref = 60, R = 1.6Q 



29813 
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FIGURE 5.1b. Inductance current, V r = 60V., R = 1.6ft 

ret 




FI GURE 5.1c. Output voltage, V re £ = 60V., R = 1.6ft 
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FIGURE S.ld. Inductance voltage for V re£ = 60V., R = 1.6fl 
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FIGURE S.3. Output voltage, V re£ = 62V. 
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FIGURE 5.5. Output voltage,, = 5 8V., 
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FIGURE 5.7. Output voltage, V r = 54V. 
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FIGURE 5.9. Output voltage, V ref s 50V. , R = 3.2ft 
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FIGURE 5.11. Output voltage, ^ re £ = 46V 



Figure 5-12 to Figure 5.1^ shows the input plot to 
the filter with average value of the pulse width approxi- 
mately corresponding to the limit cycle frequency. 

Output and input plot for the R = 1.6ft case are 
shown in Figure 5.15 to Figure 5.26. For R = 0.8ft case 
they are shown in Figure 5.27 to Figure 5*3^. 

B. ANALYSIS OF THE D-C SYSTEM 

In this paper, one of the important things to be deter- 
mined (if possible) is the exact cause of the subharmonics. 
Therefore, to get results it is necessary to make some 
change in voltage regulator system. Then it will be easy 
to see the effect of the element, (i.e., relay, thyristor). 
For this purpose the rectifier part was replaced with a 
D-C source, and the system was as shown in Figure 5-35. 

For analysis of the system the C^.5000 analog com- 
puter was used and the system output obtained for R = 0.8ft, 

R = 1.6ft, and R = 3.2ft. In this system the filter unit 
was the same as for all other tests, the only change was 
replacing the rectifier with a 65 volt battery. Fig- 
ure 5.36 shows the analog diagram for the 0^5000. Output 
and input plots for different loads are shown in Figure 
5.37 to 5.39. These outputs show that the system did not 
have subharmonics. Also these runs were made for different 
Vref conditions but the system never produced subharmonics. 
This means that if the system has a nonlinear element but 
does not have a periodic forced input then the system 
never gives subharmonics at the output. In order to have 
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FIGURE 5.13. Voltage into the filter. v re f = 56V., R 
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FIGURE 5.14. Voltage into filter, V re£ = 48V., R = 3.2ft 
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FIGURE 5.16. Output voltage, V re £ = 62V 
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FIGURE 5.18. Output voltage, V r ~ 58V 

re r 
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FIGURE 5.20. Output voltage, \ r re£ = 54V., R = 1.6S2 
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FIGURE 5.22. Output voltage, V ref - 50V., R 1 .6Q 
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FIGURE 5.24. Output voltage, = 48V. 




FIGURE 5.25. Voltage into filter, V re£ = 46V., R = 1.612 
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FIGURE 5.28. Output voltage, V r = 62V. 
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FIGURE 5.30. Output voltage, ^ TQ f = 54V. 
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FIGURE 5.33. Voltage into filter, V rg£ = 60V., R = 0 . 8« 
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FIGURE 5.34. Voltage into filter, V re£ = 50V., R = 0.8« 
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FIGURE 5.35, 



Basic block diagram of D-C case. 
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FIGURE 5.36. Analog diagram of D-C case. 
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FIGURE 5.37b. Output voltage, V re£ = 62V., R = 0 . 8« 
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FIGURE 5.38b. Output voltage, V vef = 62V., R = 1.6ft 
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FIGURE 5.39b. Output voltage, V - - 62V. 

ret 




FIGURE 5.40b. Output voltage, = 62V., R = 3.2ft 



subharmonics, a control system must satisfy two conditions, 
i.e., the system must have a periodic input forcing func- 
tion and a nonlinear element. From this it might be said 
that the exact cause of subharmonic in the limit cycling 
voltage regulator is due to the nonlinear part together 
with relay and rectifier and also the periodic input to the 
system from the AC ripple. It is important to remember 
that in order to have subharmonics the system usually must 
be lightly damped. 

The block diagram shown for D-C case at Figure 5.36's 
detail of the elements are explained in Appendix C. The 
output and input plot of system shown from Figure 5*37 to 
Figure 5 - 39 are for steady state conditions. Figure 5-^0 
shows the change in pulse duration of the relay upon the 
system reaching the steady state. 

C. GAIN CONSIDERATION 

As a last investigation of the limit cycling voltage 
regulator, the gain was changed from 1 to 10,000 and the 
output plot was obtained for different values of gain. 
Figure 5.^1, Figure 5-^2, and Figure 5*^3 show these out- 
put plots. Output shows that change in system gain did 
not effect the output wave form and it did not effect the 
subharmonic resonance. 
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FIGURE 5.41a. Output plot, V ref = 60V., R = 3.2ft, Gain - 5.0 




FIGURE 4.41b. Output plot, V re£ = 60V., R - 1.6ft, Gain = 5.0 




FIGURE 5.41c. Output plot, V re£ = 60V., R - 0.8ft, Gain 5.0 
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FIGURE 5.42a. Output plot, V re£ = 60V., R = 3.20, Gain = 1000.0 




FIGURE 5.42b. Output plot, V re£ = 60V, R = 1.60, Gain = 1000.0 




FIGURE 5.42c. Output plot, V rg£ = 60V, R = 0.80, Gain = 1000.0 



87 




FIGURE 5.43a. Output plot, V ref = 60V, R = 3.2fi, Gain = 10,000 




FIGURE 5.43b. Output plot, V ref = 60V, R = 1.6fi, Gain = 10,000 




FIGURE 5.43c. Output plot, V re£ = 60V, R = 0.8«, Gain = 10,000 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



Under some load conditions a limit cycling voltage 
regulator system gives subharmonic resonances at the out- 
put. This is not permissible in many applications. The 
goal of this thesis is to investigate the voltage regula- 
tor system and to try to find out the causes of the sub- 
harmonic resonance. 

Removal of the rectifier and substitution of a bat- 
tery removed two basic characteristics from the system: 

1. The ripple at frequency 360Hz due to the recti- 
fication of the AC supply. 

2. The delayed shut-off (i.e., continued conduction 
after the error detector signal gave the "off” signal) 
due to the thyristor firing characteristics. 

Since the subharmonics are observed to be submulti- 
ples of the limit cycle frequency, therefore, item (1) 
cannot b.e the cause of such subharmonics. Item (2) 
occurs at the limit cycle frequency, and results showed 
that the input voltage to the filter had a variable pulse 
duration for rectified AC but there was no pulse width 
modulation which probably caused the subharmonics. It 
can be seen from the input plots that the ripple begins 
at different voltage levels even though the system load 
resistance and reference voltage are held constant during 
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one analysis. The thyristor characteristics apparently 
cause this pulse width modulation. 

Another important thing is the prediction of sub- 
harmonic resonance by using the Dual input Describing 
function. If the Dual input Describing function is derived 
then by designing a compensator the subharmonics may be 
removed from the system. As a further study, the Dual 
input Describing function of the nonlinearity should be 
derived, taking into account the pulse width modulation. 
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THREE PHASE, FULL WAVE, LIMIT CYCLING VOLTAGE REGULATOR 
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APPENDIX B: DIGITAL COMPUTER PROGRAM FOR LIMIT CYCLING 

VOLTAGE REGULATOR 



The constants C(I) represent 

C(l) = Maximum value of input sine wave 
C(2) = Value of the filter inductance is henries 
C ( 3 ) = Value of the filter capacitance in farads 
C(*J) = Value of load resistance in ohms 
C ( 5 ) = Reference voltage 

C ( 6 ) = Gain of amplifier before dead zone 
C ( 7 ) = Magnitude of dead zone 
C ( 8 ) = Gain of amplifier after dead zone 
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VEL 18 



MAIN 



DATE 



70320 



19/12/10 



DIMENSION X(30) , XDOT < 30 ) , C < 1 5 ) 

C ( 10 ) =1 . 

I BP= 1 

P I =3© 141 592 
X ( 8 ) = I BP 
X< 10 )=0©0 

1 CALL INTEG2 (T,X,XDOT,C) 

X ( 1 1 ) =C ( 2 ) * XDOT ( 2 ) 

XMAX=C ( 1 ) *S IN( PI/3o ) 

THA=2©0*PI*60©0*T 

PHAS EA = ABS ( C ( 1 KSIN(THA) ) 

PHAS EB = ABS ( C ( 1 )*SI N ( THA+2 © 0*P I / 3© ) ) 

PHAS EC = ABS ( C ( 1 ) * S I N ( TH A + P I /3o ) ) 

X ( 6 ) = AMAX1 ( PHASE A, PHAS EB, PHAS EC) 
IF(PHASEAoGEoXMAX) IPH=1 
IF(PHASEB.GEoXMAX) IPH=2 
I F ( P HASECo GEo XMAX ) IPH=3 

206 E RR0R = X ( 1 ) -C ( 5 ) 

X ( 8) = I PH 
VK1=C ( 6 ) TERROR 
IF(VKloGE.C(7) ) VDZ =VK 1-C ( 7 ) 

IF(VKloLToC(7) ) VDZ=0o0 

IF( VKlcLT.CoO) VDZ=VK1 
VK2= VDZ*C ( 8 ) 

GO T0(100, 101, 102, 103, 104), IBP 

100 CONTINUE 

I F ( V K2 ) 2 07 t 207 t 20 8 

207 I BP-1 

V I N = X ( 6 ) 

X ( 5 ) =V I N 
X ( 7 ) =0o 

X(3)=(1./C(2))*(VIN-X<1)) 

X DOT ( 2 ) =X ( 3 ) 

IF(X(2)oLT o -0.01 ) X(2)=-0o01 
X(4)=(1©/C(3))*X(2)-(1©/(C(3)*C(4)))*X<1) 
XDOT ( 1 ) =X ( 4 ) 

VL AS T = V I N 
X ( 8 ) = I BP 
GO TO 1 

208 IF(VINoLToVLAST) GO TO 207 
GO T 0 ( 101,102, 103) , I PH 

101 V I N= PHASEA 
IF(VK2)207,207,105 

105 CONTINUE 
X ( 7 ) = 1 o 

X ( 5 ) = V I N 

IF(VIN-o 75) 300,301,301 

300 V I N= 0© 0 
X ( 5 ) =V I N 
I BP = 5 

GO TO 302 

301 I BP- 2 

302 X(3)=(1©/C(2))*<VIN-X<1)) 

IF(X (2)®LT©-0©01 ) X ( 3 ) = 0© 0 
X DOT ( 2 ) =X ( 3 ) 

. IF(X(2)©LT©-0©01 ) X ( 2 ) =-0© 01 

X<4)=(1©/C(3))*X(2)-<1©/(C<3)*C<4)>)*X<1) 
XDOT ( 1 ) =X ( 4 ) 

VLAST=VI N 
X ( 8 ) = I BP 
GO TO 1 

102 VIN=PHASEB 

I F(VK2 ) 207, 207 , 106 

106 CONTINUE 
X ( 7 ) =2 © 

X ( 5 ) =V I N 

IF(VIN-© 75) 3 03,304,304 

303 V I N= 0© 0 
X ( 5 ) =V I N 
I BP= 5 

GO TO 305 

304 I BP= 3 

305 X(3)=(1©/C(2))*(VIN-X(1)) ■ 
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.EVEl 18 



MAIN 



DATE 



70320 



19/12/10 



I F ( X ( 2 ULTo-OoOl ) X ( 3 ) = 0o 0 
XDOT ( 2 ) =X ( 3 1 

IF(X(2)«LTo-0.01 ) X(2)=-0o01 
X(4)=(lo/C(3))*X(2)“(l©/(C(3)*C(4) ) ) *X ( 1 ) 
XDOT ( 1 ) =X ( 4 ) 

VLAST=VI N 
X ( 8 ) = I BP 
GO TO 1 

103 V IN=PHASEC 
IF(VK2)207,207,107 

107 CONTINUE 
X ( 7 ) =3® 

X ( 5 ) =V I N 

IF(VIN-„ 75)306,307,307 

306 V I N= Oo 0 
X ( 5 ) =V I N 
I 8P = 5 

GO TO 308 

307 I BP=4 

308 X ( 3 ) =(lo/C(2) )*(VIN-X(1) ) 

IF(X (2)oLTo~0*01 ) X( 3)=0o0 
X DOT ( 2 ) =X ( 3 ) 

IF(X(2)oLTo-*0®01 ) X(2)=-0o01 
X(4)=(lo/C(3) ) *X ( 2 )-(lo/(C(3)**C(4) ) ) *X ( 1 ) 
XDOT ( 1 ) =X ( 4 ) 

VLAST=V I N 
X ( 8 ) = I BP 
GO TO 1 

104 V IN = Oo 0 
X ( 7 ) =4o 

X ( 5 ) =V I N 

X(3)=(l*/C(2) ) * ( V I N- X ( 1) ) 

IF(X(2).LTo-0«01 ) X ( 3 ) = Oo 0 
XDOT ( 2 ) =X ( 3 ) 

IF(X(2) o LTo-0.01 ) X(2)=-0o01 
X(4)=(l./C(3))*X(2)-(lo/(C(3)*C(4)))*X(l) 
XDOT ( 1 ) =X ( 4 ) 

I BP= 5 

IF(VK2 1309,309,310 

309 I BP= 1 

310 CONTINUE 
X ( 8 ) = I BP 
GO TO 1 
END 



94 



GULER M.Ao FORCED LIMIT CYCLING REGULATOR 
ONE RUN IS CALLED FOR 



N04 



INPUT DATA RECORD 
ORDER OF EQUATIONS = 2 



THE NON-ZERO CONSTANTS, C(I), ARE 



C < I ) = Oo 65C0E 02 
C< 2) = Oo 5000E-03 
C( 3) = OolOOOE 00 
C< 4) = O 0 I 6 OOE 01 
C( 5) = O 06 OOOE 02 
C ( 6 ) = OolOOOE 01 
C< 8) = OolOOOE 01 



ALL THE INITIAL CONDITIONS ARE ZERO 

THE COLUMN HEADINGS AND THE CORRESPONDING VARIABLES ARE 



NO GRAPHS ARE REQUIRED 



INITIAL TIME 
FINAL TIME 
STEP SIZE 



OoO 

OolOOOE 00 
0 o 22 22 E-04 



VIN 

VOUT 

CURRT 
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GULER Mo A® FORCED LIMIT CYCLING REGULATOR NG4 



TIME VIN 



VOUT CURRT 



Oo 0 


Oo 562 92E 


02 


OoO 




Oo 44440E-03 


0 o 60 923 E 


02 


0oll435E 


00 


Oo 88880E-03 


0 o 63848E 


02 


Oo 467 59E 


00 


Oo 13332E-02 


Oo 64 986E 


02 


Oo 106 92 E 


01 


Oo 17776E-02 


Oo 64303E 


02 


Oo 192 05E 


01 


Oo 2 2220E- 02 


0 « 6 1 8 20E 


02 


Oo 301 49 E 


01 


Oo 26664E-02 


Oo 576C6E 


02 


0 o 43375E 


01 


Oo 311 08E-02 


Oo 59918E 


02 


Oo 5 86 89 E 


01 


0o35552E-02 


0 o 63 28 1 E 


02 


Oo 760 8 8 E 


01 


Oo 39996E-02 


0 o 64 871 E 


02 


Oo 95628 E 


01 


Oo 44439 E-02 


Oo 64 645 E 


02 


0 0 1 17 29E 


02 


Oo 48883E-02 


Oo 6261 OE 


02 


Oo 140 97 E 


02 


Oo 53326E- 02 


Oo 588 21E 


02 


Oo 1 6649E 


02 


Oo 57770E-02 


Oo 588 C5E 


02 


Oo 193 61 E 


02 


Oo 6221 3 E-02 


Oo 62 6 COE 


02 


0 0 222 2 1 E 


02 


Oo 66657E-02 


0 o 64 64 IE 


02 


Oo 2 52 3 1 E 


02 


Oo 7 1 10 1 E-02 


0.64873E 


02 


0o28388E 


02 


Oo 75544E-02 


0 o 632 8 9E 


02 


Oo 31680E 


02 


Oo 79988 E-02 


0 o 59933E 


02 


Oo 350 86E 


02 


Oo 8443 IE-02 


Oo 57588E 


02 


Oo 385 8 1 E 


02 


Oo 88875E-02 


O06I8C8E 


02 


0 o42 1 43 E 


02 


Oo 9331 8 E-02 


0 0 64 2 ^ 7 E 


02 


0 0 457 73 E 


02 


Oo 97762E-02 


Oo 64987E 


02 


Oo 4°4 64E 


02 


Oo 10221E-01 


0 0 63 85 7E 


02 


0 0 532 07E 


02 


Oc 10665E-01 


Oo 60 939E 


02 


Oo 56980E 


02 


Oo 11 109E-01 


Oo 5631 5E 


02 


Oo 607 57E 


02 


Oo 11554E-01 


Oo 50 1 1 5E 


02 


0 0 645 07E 


02 


Oo 11998E-01 


Oo 42 5 1 IE 


02 


O068I88E 


02 


Oo 1 2442E-01 


0o33717E 


02 


0 o 717 58E 


02 


0oi2887E-01 


0 0 239 80E 


02 


0 o 75167E 


02 


Oo 13331E-01 


0ol3571E 


02 


Oo 783 65 E 


02 


0ol3775E-01 


0 0 27 820E 


01 


0 o 8 12 99 E 


02 


Oo 14220E-01 


OoO 




Oo 8391 7 E 


02 


Oo 14664E-01 


OoO 




Oo 86198E 


02 


Oo 15108E-01 


OoO 




0o88132E 


02 


Oo 1 5 553E-01 


OoO 




Oo 897 1 4E 


02 


Oo 1 5997 E- 01 


OoO 




Oo 909 37 E 


02 


Oo 16441E-01 


OoO 




Oo 917 98E 


02 


Oo 1 6886E- 01 


OoO 




Oo 922 95 E 


02 


Oo 17330E-01 


OoO 




Oo 924 27 E 


02 


Oo 1 7775E-01 


OoO 




Oo 922 28 E 


02 


Oo 1 82 19E-01 


OoO 




Oo 9 19 72 E 


02 


Oo 18663E— 01 


OoO 




Oo 917 17 E 


02 


Oo 19108E-01 


OoO 




Oo 91462E 


02 


Oo 19552E-01 


OoO 




Oo 9 12 08 E 


02 


Oo 19996E— 01 


OoO 




Oo 909 55 E 


02 


Oo 2044 1 E-0 1 


OoO 




Oo 907 03E 


02 


Oo 2088 5E-01 


OoO 




Oo 904 5 1 E 


02 


Oo 2 1329E-01 


OoO 




Oo 902 00E 


02 


Oo 2 1774E-01 


OoO 




Oo 899 49 E 


02 



OoO 

Oo 52 178E 02 
0 © 107 5 2E 03 
Oo 16424E 03 
Oo 2205 2E 03 
Oo 2745 2E 03 
Oo 32447E 03 
Oo 37 1 3 8E 03 
Oo 42028E 03 
Oo 46 97 5E 03 

Oo 5 1 799E 03 
Oo 5632 1 E 03 
Oo 60365E 03 
Oo 63884E 03 
Oo 6744 5E 03 
Oo 71 0 05E 03 
Co 74392E 03 
Oo 77 432E 03 
Oo 79954E 03 
Oo 8 1 82 6E 03 

0o83557E 03 
Oo 85267E 03 
Oo 867 94E 03 
Oo 8797 IE 03 
Go 88633E 03 
Oo 88623E 03 
Oo 87797E 03 
Oo 8602 5E 03 
Oo 832 00E 03 
Oo 79239E 03 

Oo 74087E 03 
Oo 67717E 03 
Oo 60 4 02E 03 
Oo 52839E 03 
Oo 45 08 9E 03 
Oo 37 1 8 3E 03 
Oo 29 1 5 2E 03 
Oo 21 028E 03 
Oo 12 845 E 03 
Oo 46328E 02 

- Oo IOC OOE-0 1 
-Oo 10000E-01 
-Oo 100 00E-01 
-Oo 100 00E-01 
-Oo 10 OOOE-0 1 
-Oo 10000E-01 
-Oo 10000E-01 
-Oo 100 OOE-0 1 
-Oo 100 00E-01 
-Oo 100 00E-01 
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GULER Mo Ao FORCED LIMIT CYCLING REGULATOR 



TIME 


VIN 


VOUT 


CURRT 


Oo 2221 8E-01 


OoO 


Oo 897 00E 


02 -Oo 100 OOE-0 1 


0o22662E-01 


OoO 


Oo 89451 E 


02 -Oo 10 0 00E-01 


Oo 23 1 07E-01 


OoO 


Oo 892 02E 


02 -Oo 10000E-01 


Oo 23551 E-01 


OoO 


0 o 8 89 55 E 


02 -Oo 100 00E-01 


Oo 23995E-01 


OoO 


Oo 88708E 


02 -Oo 100 00 E-0 1 


Oo 24440E-01 


OoO 


Oo 88462 E 


02 -Oo 100 00E-01 


Oo 24884E-01 


OoO 


0 o 8 82 1 6E 


02 -Oo 100 OOE- 0 1 


Oo 25329E-01 


OoO 


0 o 879 7 1 E 


02 -Oo 10000E-01 


Oo 25773E-01 


OoO 


Oo 877 27E 


02 -Oo 100 00E-01 


Oo 262 1 7 E-01 


OoO 


0 o 87484E 


02 - Oo 1 0000E-01 


Oo 26662E-01 


OoO 


Oo 87241 E 


02 -Oo 1 0000E-01 


Oo 27 1 06 E- 01 


OoO 


Oo 869 99 E 


02 -Oo 10000E-01 


Oo 27550E-01 


OoO 


Oo 86757E 


02 -Oo 10000E-01 


Oo 27995E-01 


OoO 


Oo 865 1 6E 


02 -Oo 100 00E-01 


Oo 2 8439 E- 01 


OoO 


0 o 862 76E 


02 -Oo 10000E-01 


Oo 28883E-01 


OoO 


Oo 8603 7 E 


02 - Oo 10 0 OOE-0 1 


Oo 2932 BE- 01 


OoO 


Oo 857 98 E 


02 -Oo 100 OOE-0 1 


Oo 297 72 E-01 


OoO 


Oo 855 60 E 


02 -Oo 10000E-01 


Oo 30216E-01 


OoO 


0.85322E 


02 - Oo 100 00E-01 


Oo 3 066 1 E-0 1 


OoO 


Oo 8 50 8 5 E 


02 - Oo 1 00 OOE-0 1 


0o31105E-01 


OcO 


Oo 848 49 E 


02 -Oo 10000E-01 


Oo 3 1 549E-01 


OcO 


0 © 8 46 1 3 E 


02 -Oo 10000E-01 


Oo 3 1 994E-01 


OoO 


Oo 843 79 E 


02 -Oo 10000E-01 


Oo 3243 8 E- 01 


OoO 


Oo 84 1 44E 


02 - Oo 10 0 00 E -01 


Oo 328S3E-01 


OoO 


0 o 839 1 1 E 


02 -Oo 1 00 00 E-01 


Oo 33327E-01 


OoO 


Oo 83678E 


02 -Oo 10000E-01 


Oo 33771 E-01 


OoO 


0.83446E 


02 -Oo 1 00 OOE -0 1 


Oo 342 1 6E-01 


OoO 


Oo 832 1 4E 


02 -Oo 10000E-01 


Oo 34660 E-01 


OoO 


Oo 829 83 E 


02 -Oo 100 00E-01 


Oo 35 104E-01 


OoO 


Oo 827 53E 


02 -Oo 10000E-01 


Oo 35549E-01 


OoO 


Oo 82 523 E 


02 -Oo 10000E-01 


Oo 35993E-01 


OoO 


Oo 822 94E 


02 -Oo 1 0000 E-01 


Oo 36437E-01 


OoO 


Oo 820 65 E 


02 -Oo 10000E-0 1 


Oo 3 6882 E- 01 


OoO 


Oo 81 837E 


02 -Oo 10000E-01 


Oo 37326E-01 


OoO 


O 08 I 6 IOE 


02 -Oo 10000E-01 


Oo 37770E-01 


OoO 


0o81384E 


02 -Oo 100 00E-01 


Oo 382 1 5E-0 1 


0.0 


Oo 81 1 58E 


02 -Oo 100 00E-01 


Oo 38659E-01 


OoO 


Oo 809 32 E 


02 -Oo 10000E-01 


Oo 39103E-01 


0.0 0 


Oo 80708E 


02 -Oo 100 00E-01 


Oo 39548E-01 


OoO 


Oo 80484E 


02 -Oo 10000E-01 


Oo 39992 E-01 


OoO 


0o80260E 


02 -Oo 10 0 00E-01 


0 o 40436E-01 


OoO 


0® 80038E 


02 -Oo 10000E-01 


Oo 4088 1 E-01 


OoO 


Oo 798 1 5 E 


02 -Oo 100 OOE- 0 1 


Oo 41 32 5E-01 


OoO 


Oo 795 94 E 


02 -Oo 100 00E-01 


Oo 41770E-01 


OoO 


Oo 793 73E 


02 -Oo 10000E-01 


Oo 422 14E-01 


OoO 


Oo 79 1 52E 


02 -Oo 10000E-01 


Oo 42658E-01 


OoO 


Oo 789 33E 


02 -OolOOOOE-Ol 


Oo 43 1 03E-01 


OoO 


Oo 7 87 1 4E 


02 -Oo 10000E-01 


Oo 43 547E-01 


OoO 


Oo 784 95 E 


02 -Oo 10000E-01 


Oo 4399 IE-01 


OoO 


Oo 782 77 E 


02 -Oo 10000E-01 



N04 
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FORCED LIMIT CYCLING REGULATOR 



TIME 



Oo 44436E- 01 
Oo 44880E-01 
Oo 45 324E-01 
Oo 45769 Ef- 01 
0o46213E-01 
Oo 46657E-01 
Oo 47 1 02 E- 01 
Oo 47546E-01 
Oo 47990E-01 
Oo 48435E-01 

Oo 4 8879E- 01 
Oo 49324E- 01 
Oo 4 9768 E-0 1 
Oo 502 1 2 E-01 
Oo 50657E-01 
Oo 51101E-01 
Oo 51 545E-01 
Oo 5 1 990E-0 1 
Oo 52434E-01 
0o52878E-01 

Oo 53323E-01 
0© 537 67 E-01 
Oo 5421 IE-01 
Oo 54656E-01 
0 o 55100E-01 
Oo 55544E-01 
Oo 55989E-01 
Oo 5 643 3E-01 
Oo 56878E-01 
Oo 57 322 E-01 

Oo 57766E-01 
Oo 582 1 1 E-01 
Oo 5 8655E-01 
Oo 5 9099E-01 
Oo 59544E-01 
Oo 5 9988E-01 
Oo 60432E-01 
Oo 60877 E- 01 
Oo 6 132 1 E-01 
Oo 61765E-01 

Oo 622 10 E-01 

Oo 62 654E-01 
Oo 63097 E-01 
Oo 63541 E-01 
Oo 63984E- 01 
Oo 6442 8E-01 
Oo 64871 E-01 
0o65315E-01 
Oo 65758E-01 
Oo 662 02 E-01 



GULER Mo Ao 



VOUT 




Oo 78060E 


02 


0.77843E 


02 


Oo 77627 E 


02 


Oo 774 12F 


02 


0.77197E 


02 


0 o 769 82 E 


02 


0.76769E 


02 


Oo 765 55 E 


02 


Oo 76343E 


02 


0o76131E 


02 


Oo 75920E 


02 


Oo 757 09E 


02 


Oo 75499E 


02 


Oo 7 52 89E 


02 


Oo 750 80E 


02 


Oo 748 72 E 


02 


Oo 746 64E 


02 


0«74456E 


02 


Oo 742 5 OE 


02 


Oo 74044E 


02 


Oo 738 38E 


02 


Oo 73633E 


02 


Oo 73429E 


02 


Oo 732 2 5E 


02 


Oo 7302 1 E 


02 


Oo 728 1 9E 


02 


0.72617E 


02 


0 o 724 1 5E 


02 


Oo 722 14E 


02 


Oo 720 13E 


02 


Oo 718 13E 


02 


Oo 716 14E 


02 


Oo 714 15E 


02 


Oo 712 17E 


02 


0 o 71019E 


02 


Oo 70822E 


02 


Oo 70625 E 


02 


Oo 70429E 


02 


Oo 702 34E 


02 


Oo 70039E 


02 


Oo 698 44 E 


02 


Oo 696 5 OE 


02 


0 o 694 57 E 


02 


0 o 6 92 64E 


02 


Oo 690 72 E 


02 


O06888OE 


02 


Oo 68689E 


02 


0 0 6 84 9 8E 


02 


0 0 6 83 0 8 E 


02 


O068II8E 


02 



CURRT 



-Oo 10000E-01 
-OolOOOOE-Ol 
-Oo 100 00 E-0 1 
-Oo 10000E-01 
-Oo 10 0 OOE-0 1 
-Oo 10000E-01 
-Oo 100 00E-01 
-Oo 100 OOE-0 1 
-Oo 100 00E-01 
-Oo 100 OOE-0 1 

-Oo 100 00E-01 
-Oo 10000E-01 
-Oo 10000 E-0 1 
-Oo 10000E-01 
-Oo 100 00 E-01 
-Oo 10000E-01 
-Oo 10000E-01 
-Oo 100 OOE-0 1 
-Oo 10000E-01 
-Oo 10000E-01 

-Oo 10000E-01 
-Oo 10000E-01 
-Oo 10000E-01 
-Oo 10000E-01 
-Oo 1 00 OOE-0 1 
-Oo 10000E-01 
-Oo 100 00E-01 
-Oo 10 0 OOE-0 1 
-Oo 100 OOE-0 1 
-Oo 10000E-01 

-Oo 100 00E-01 
-Oo 100 00E-01 
-Oo 100 OOE-0 1 
-Oo 10000E-01 
-Oo 100 OOE-O 1 
-Oo 100 OOE-O 1 
-Oo 10000E-01 
-Oo 10000E-01 
-Oo 10000E-01 
-Oo 10000E-01 

-Oo 10000E-01 
-OolOOOOE-Ol 
-Oo 10000E-01 
-Oo 100 00E-01 
-Oo 10000E-01 
-Oo 100 00 E-01 
-Oo 10000E-01 
-Oo 100 00E-01 
-Oo 10 0 OOE-O 1 
-Oo 10000E-01 



VIN 



OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

0 c 0 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 

OoO 



N04 
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GULER Mo Ao FORCED LIMIT CYCLING REGULATOR 



TIME 


VIN 




VOUT 




CURRT 


Oo 66645E-01 


OoO 




Oo 67929E 


02 


-Oo 10000E-01 


Oo 67 08 8 E- 01 


OoO 




0. 6774 1 E 


02 


-Oo 100 00E-01 


Oo 67 532E-01 


OoO 




0« 675 52E 


02 


-Oo 10000E-01 


Oo 679 75 E- 01 


OoO 




Oo 673 65E 


02 


-Oo 100 00E-01 


Oo 68419E-01 


OoO 




0.67178E 


02 


-Oo 100 00E-01 


0 o 68862E-01 


OoO 




Oo 669 91 E 


02 


-Oo 10000E-01 


Oo 69306E-01 


OoO 




Oo 668 05 E 


02 


-Oo 10000E-01 


Oo 69749E-01 


OoO 




0« 666 2 OE 


02 


-Oo 1 00 00 E-0 1 


Oo 70 1 93 E-01 


OoO 




Oo 66435 E 


02 


-Oo 10000E-01 


Oo 7 063 6 E- 01 


0.0 




0. 662 50 E 


02 


-Oo 100 00 E -01 


Oo 71080E-01 


OoO 




Oo 660 66 E 


02 


-Oo 10000E-01 


Oo 71 523E-01 


OoO 




Oo 658 83E 


02 


-Oo 10000E-01 


Oo 7 1966 E-0 1 


OoO 




0 . 6 57 00 E 


02 


-Oo 10000E-01 


Oo 7241 0E-01 


OoO 




Oo 655 1 8E 


02 


-Oo 10000E-01 


Oo 72853E-01 


OoO 




Oo 653 3 6 E 


02 


-Oo 100 OOE-Ol 


Oo 73297 E-01 


OoO 




Oo 6 5 1 54E 


02 


-Oo 10000E-01 


Oo 73 740E- 01 


OoO 




Oo 649 73 E 


02 


-Oo 10000E-01 


Oo 74 1 84E- 01 


OoO 




0.64793E 


02 


-Oo 100 00E-01 


Oo 74627E-01 


OoO 




Oo 646 1 3 E 


02 


-Oo 10000E-01 


Oo 75071 E-01 


OoO 




Oo 64434E 


02 


-Oo 100 00E-01 


Oo 75 51 4E-01 


OoO 




Oo 642 55 E 


02 


-Oo 100 00 E-0 1 


Oo 75958E-01 


OoO 




0. 64076E 


02 


-Oo 10000E-01 


Oo 7 6401 E-01 


OoO 




0.63898E 


02 


-Oo 10000E-01 


Oo 76845 E-01 


OoO 




Oo 637 2 IE 


02 


-Oo 10000E-0 1 


Oo 77288E-01 


OoO 




Oo 63544E 


02 


-Oo 10 0 00 E-01 


Oo 7773 1 E-01 


OoO 




0 o 6 33 6 8 E 


02 


-Go 10000E-01 


Oo 7 8 1 75E-01 


OoO 




0© 63 1 92 E 


02 


-Oo 10 0 00E-01 


Oo 78618E-01 


OcO 




0© 63 0 1 6E 


02 


-Oo 1 00 00E-01 


Oo 7 90t>2 E-0 1 


OoO 




Oo 628 41 E 


02 


-Oo 10 0 00 E-0 1 


Oo 79505E — 01 


OoO 




Oo 62667E 


02 


-Oo 10000E-01 


Oo 79949E-01 


OoO 




Oo 624 93 E 


02 


-Oo 10000E-01 


Oo 80392 E-01 


OoO 




Oo 623 1 9 E 


02 


-Oo 10000E-01 


Oo 80836 E- 01 


OoO 




Oo 621 46E 


02 


-0© 10000E-01 


Oo 81279E-01 


OoO 




0o61974E 


02 


-Oo 10000E-01 


Oo 8 1 72 3 E- 01 


OoO 




0. 61 8 02 E 


02 


-Oo 10000E-01 


Oo 8 2 1 66 E-01 


OoO 




Oo 6 16 30E 


02 


-Oo 10000E-01 


Oo 82609E-01 


OoO 




0.61459E 


02 


-Oo 10000E-0 1 


Oo 83053E-01 


OoO 




0«61288E 


02 


-Oo 100 00E-01 


Oo 8 3496 E-01 


OoO 




O 06 III 8 E 


02 


-Oo 10000E-01 


Oo 83 940E- 01 


OoO 




Oo 60948E 


02 


-Oo 10 0 00 E-01 


Oo 84383E-01 


OoO 




Oo 607 79E 


02 


-Oo IOC OOE-Ol 


Oo 84827E-01 


OoO 




O 06 O 6 IOE 


02 


-Oo 100 OOE-Ol 


Oo 85270E-01 


OoO 




Oo 60442E 


02 


-Oo 10000E-01 


Oo 857 14E-01 


OoO 




Oo 602 74E 


02 


-Oo 10000E-01 


Oo 861 57E-01 


OoO 




0.60107E 


02 


-Oo 100 OOE-Ol 


Oo 866 0 1 E-0 1 


0o61298E 


02 


Oo 59940E 


02 


-Oo 10000E-01 


Oo 87044E-01 


Oo 64042E 


02 


Oo 59776E 


02 


Oo 1 9466E 01 


Oo 8748 8E-01 


Oo 64999E 


02 


0 0 596 2 8 E 


02 


Oo 63 648E 01 


Oo 8793 1 E-0 1 


0 o 64 144E 


02 


Oo 59501 E 


02 


0ol0951E 02 


Oo 88374E-01 


Oo 61 501E 


02 


Oo 593 92 E 


02 


Oo 14083E 02 



N04 
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FORCED LIMIT CYCLING REGULATOR N04 



GULER MoA. 



TIME 


VIN 




VOUT 




CURRT 




O0888I8E-OI 


0 0 57 14 2E 


02 


Oo 592 91 E 


02 


Oo 14188E 


02 


Oo 89261E-01 


Oo 602 8 5E 


02 


Oo 59186E 


02 


Oo 13 1 91 E 


02 


Oo 897 05E-01 


0.63489E 


02 


0® 59085E 


02 


0® 15766E 


02 


Oo 901 48 E-0 1 


Oo 6492 2E 


02 


Oo 590 01 E 


02 


Oo 20489E 


02 


Oo 90592 E- 0 1 


Oo 64546E 


02 


0® 58940E 


02 


3o25748E 


02 


0o91035E~01 


Oo 62 369E 


02 


Oo 58901E 


02 


Oo 2 991 4E 


02 


Oo 91479E-01 


Oo 58454E 


02 


Oo 5 88 74E 


02 


Oo 31394E 


02 


Oo 9 1922E-01 


0 0 59 1 5 5E 


02 


Oo 58847E 


02 


Oo 30 28 1 E 


02 


Oo 92366E-01 


0 « 62 8 1 3 E 


02 


0® 58822E 


02 


Oo 32318E 


02 


Oo 92809E- 01 


Oo 6471 9E 


02 


Oo 5 88 1 1 E 


02 


Oo 3 685 IE 


02 


Oo 93252 E-01 


0 0 64 82 1 E 


02 


0® 588 23E 


02 


Oo 42278E 


02 


Oo 93696E-01 


0 « 63 1 1 6E 


02 


Oo 5 88 5 9E 


02 


Oo 46969E 


02 


Oo 941 39E-01 


0 0 59 65 OE 


02 


0 0 5 89 1 OE 


02 


Oo 493 1 8E 


02 


Oo 94583 E-01 


Oo 57 909E 


02 


Oo 58964E 


02 


Oo 48 260 E 


02 


Oo 95026E-01 


Oo 62 01 4E 


02 


Oo 59015E 


02 


Oo 49249E 


02 


Oo 95470E- 01 


0 0 64 3 90E 


02 


0® 59078E 


02 


Oo 53074E 


02 


Oo 9591 3E-01 


Oo 649 71 E 


02 


Oo 591 60E 


02 


Oo 58153E 


02 


Oo 96357E- 01 


Oo 63739E 


02 


Oo 59265E 


02 


Oo 6285 8E 


02 


Oo 96800E-01 


0o60731E 


02 


Oo 59387E 


02 


Oo 65 574E 


02 


0o97244E-01 


0 0 56550E 


02 


Oo 59513E 


02 


Oo 64756E 


02 


Oo 97687 E-01 


0 0 61 094E 


02 


0.59632E 


02 


0o64211E 


02 


Oo 98131 E-01 


Oo 63 936F 


02 


Oo 597 56E 


02 


Oo 66849E 


02 


Oo 98574E-01 


Oo 64994E 


02 


0 0 598 96E 


02 


Oo 71 1 06E 


02 


Oo 9901 7 E-01 


Oo 64 2 3 9 E 


02 


0o60055E 


02 


Oo 75 366 E 


02 


Oo 99461 E-01 


Oo 61 694E 


02 


0o60230E 


02 


Oo 78007E 


02 



AT 4500 INTEGRATION STEPS 
ONE RUN CALLED FOR HAS BEEN COMPLETED., 
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APPENDIX C 



SUMMER AMPLIFIER: 

Amplifiers on the patch board are addressed as an 
even number A000 to A12 6 , and are represented with a gain 
symbol at the input to the amplifier 




INVERTERS : 

Sign changes are addressed from A200 to A217 and are 
represented as 



INPUT 




OUTPUT 



INTEGRATOR : 

These amplifiers are addressed on the patch board as 
odd numbers (A001 to A0^7) and are integrators. They are 
represented with a symbol 




A007 and A011 were used. 
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COMPARATOR : 



A kind of device which gives Ov . or -6v. logic level 
according to the comparison of the inputs. Represented 
with symbol 



X 


Cot 


Cot 


Y 


1 


CoT 

> 



such that if 

( C = -6v. 

X+Y >o 1 — 

/ C = Ov. 

L a 

f C = Ov. 

\ a 

X+Y<o \ — 

/ C -6v. 

V a 

DIGITAL/ANALOG SWITCH (DA): 

It is a switch and is represented with a symbol 
| LOGIC 




The switch is closed if logic signal is -6v. and open if 
logic signal is Ov. These logic levels usually come from 
a comparator. Switching time is approximately equal to 
1 ysec. 
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